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FOREWORD 

This report represents the results of work performed by 
the Thermal Environment Section of the Aero -Mechanics Depart- 
ment of the Lockheed Missiles & Space Company, Huntsville 
Research & Engineering Center, for the National Aeronautics and 
Space Administration, Marshall Space Flight Center, Huntsville, 
Alabama under contract NAS8-21347. The NASA contract moni- 
tor is John G. Austin, Jr ., of the MSFC Astronautics Laboratory. 
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SUMMARY 

Two cryogenic storage tanks 20 and 105 inches in diameter, respec- 
tively, were analyzed using a Lockheed-developed Thermal Analyzer Com- 
puter program. This program is suited to the thermal analysis of venting 
cryogenic tanks in that it incorporates not only the natural convection along 
the walls of the tanks but also the thermal effects of transfer of the venting 
gas. For the purpose of comparison, three methods were used in analyzing 
each of the tanks. These methods were: 

• The venting gas was assumed to be stationary, and its temper- 
ature was constrained only at the liquid surface and at the 

top of the vent line. 

• The venting gas was assumed to be stationary, and its temper- 
■ ature was fixed at the evaporation temperature at all points. 

• The venting gas was assumed to be flowing out the vent line at 
the rate dictated by boiloff, and its temperature was constrained, 
only at the liquid surface and the top of the vent line. 

For the 20-inch tank the heat leak to the liquid for the three methods of anal- 
ysis were 7.58, 4. 15 and 5.79 3tu/hr, respectively. Test data for the 20-inch 
tank revealed a mass boiloff rate of 0.0508 lbm/hr which corresponds to a 
heat leak to the liquid of 4.38 Btu/hr. 

For the 105-inch tank, the heat leak to the liquid for the three methods 
of analysis were 133.37, 79.52 and 85.99 Btu/hr, respectively. No test data 
are yet available. 

These results indicate that the approximations made in the analysis con- 
cerning the venting gas has a significant effect on the predicted heat rate. The 
capability of treating the thermal effects of mass flow is a needed tool in anal- 
yzing venting cryogen tanks. 
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Section 1 
INTRODUCTION 

The true effectiveness of high performance insulation (HPI) in actual 
applications cannot be determined simply' by idealized calorimeter tests of 
the materials. In actual applications the thermal conductivity tends to be 
greater than the predicted value because of gaps, joints, penetrations, un- 
controllable compression, offgassing and other effects. Therefore, full 
and subscale tests of insulated cryogenic tanks are being conducted. How- 
ever, costs would be excessive to conduct tests for all future cryogenic tank 
applications. Therefore, it will be advantageous if adequate analytical, 
rather than experimental, prediction techniques could be used. Before an 
analytical prediction could be considered adequate, however, it is necessary 
to conduct comparisons between analytical and experimental predictions of 
the same tanks. If the analytical predictions do indeed match test data for a 
variety of tank designs, then, in the future, confidence could be placed in the 
predictions alone, making them a more valuable tool for use in design work. 
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Section 2 

TECHNICAL DISCUSSION 

2.1 THE NASA/MSFC 20-INCH CRYOGENIC TANK 

2.1.1 Experimental Results 

The NASA/MSFC 20-in. cryogenic tank, shown schematically, in Fig.l, 
was tested in a vacuum facility at MSFC. The test was documented in Ref. 1, 
and the following information is found therein. The tank was insulated with 
an average of 1,37 in. of l/ 4- mil double-aluminized Mylar and red polyure- 
thane foam insulation. After 96 hours of testing, the nitrogen gas boiloff was 

steady at 0.0508 Ibm/hr. The chamber pressure was estimated to be between 
-7-5 

1x10 and 1x10 torr Hg. The insulation backside pressure was esti- 

-3 -3 

mated to be between 1 x 10 and 50 x 10 torr Hg. The boiloff gas was 
maintained at slightly above atmospheric pressure 15.25 psia) by a throttle 
valve and was measured to be room temperature at the vent exit. 

Based on the steady state boiloff rate the heat rate to the nitrogen can 
be determined. 


Q = mh fg = 0.0508 x 86.2 = 4. 38 Btu/hr. 


2.1.2 Analytical Results 

The NASA/MSFC 20-inch cryogenic tank was analyzed using a Lockheed- 
developed Thermal Analyzer computer program. Figure 1 shows a schematic 
of the tank. The liquid level was chosen by setting the gas fraction in the top 
of the tank at approximately 5% which is consistent with the experimental 
data. Only the portions of the tank above the adiabatic separator, including 


2 
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the neck, were analyzed. The heat flow in regions of the tank below the adia- 
batic separator was assumed to be one -dimensional and, therefore, not in 
need of computer analysis. Since there are no penetrations to the tank other 
than the neck, heat flow in the portions of the tank below the separator are 
assumed to be perpendicular to the insulation layers and the tank wall itself. 
No heat is expected to flow down the tank wall since the uniform liquid tem- 
perature holds the thin-walled tank at a uniform temperature. Therefore, 
the heat rate to the liquid below the separator would be 


q = k A 


AT 

AX 


where 


k = the insulation's effective thermal conductivity 

A = the insulation area (below the separator) 

AT = the temperature differential through the 
insulation 

AX = the insulation thickness 


The critical item in attaining the correct q through that portion of the insu- 
lation is that the proper value for k be chosen. The other parameters A, 

AT and AX are known within close tolerence relative to k and are listed 
below: 

2 

A = 15.41 ft (calculated from tank dimensions) 

AT = 70° - (- 320) = 390°F 
AX = 1.37 in. 

The insulation used on the tank test was the l/4-mil double aluminized 
Mylar and red polyurethane foam material. At the average temperature of 
the insulation (-12 5°F) the experimental value of thermal conductivity is 
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5 _+ 1 x 10 ^ Btu/hr-ft-°F. This results in a heat flux to the liquid below the 
separator of 2.19 rate jf0.44 Btu/hr. 

The remaining heat rate to the liquid is that entering above the 
separator which includes flow down the tank wall, down the fill line and con- 
duction down the gas. All these heat rates are calculated in the thermal 
analysis. The analysis was performed assuming that the tank was insulated 
with 1.37 in. of the double aluminized Mylar and red polyurethane foam HPI. 
The thermal conductivity normal to the insulation was assumed to be 5 4 1 x 10 
Btu/hr-ft-°F while the thermal conductivity parallel to the insulation was 
assumed to be 1.55x10 ^ Btu/hr -ft- °F . Temperature-dependent thermal con- 
ductivity values for the aluminum tank, the stainless steel fill and vent lines 
and the nitrogen gas were used in the analysis. All the heat flow paths in the 
analysis were conduction except those convection heat paths connecting the gas 
to the tank wall and to the fill and veht lines . The convection heat flow paths 
(thermal resistors) were 

R = i (2) 

conv hA 

where 

h = the convection heat transfer coefficient 

A = heat flow area 

The heat transfer coefficient was obtained from the laminar and turbu- 
lent natural convection equations on a vertical wall. (Ref. 2) 

Laminar 

/ Pr \l/4 

N “ X = °- 508 (o^sgfp? Gr x Pr ) < 3 > 
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Turbulent 


Nu 

x 


0.0295 (Gr x ) 2 ^ 5 Pr 7 ^ 15 

275 




+ 0.495 Pr 


2/3 y 


( 4 ) 


In the test, the tank contained liquid nitrogen which boils at - 320°F . The 
gas properties were evaluated for gaseous nitrogen at -280°F. A value some- 
what higher than -320°F was chosen since the gas temperature will increase as 
it passes up the vent lines* The error caused by not using temperature depen- 
dent gas properties is insignificant. Therefore, the temperature at which the 
properties are evaluated is not critical. The resulting equations were: 


Gr = 3.45 x 10 8 (T + 320) x 3 
x ' w ' 


(5) 


Laminar 

h = 0.294(T +320) 1 / 4 x" 1 ^ 4 Btu/hr-ft 2 -°F (6) 

w 


Turbulent 


h = 0.323(T +320) 2 / 5 x 1 / 5 Btu/hr-ft 2 -°F 
w 


(7) 


q 

The transition Gr number was set at 10 in the analysis. In both of these equa- 
tions, h is dependent on the surface temperature, T, and the distance, x, along 
the surface from the leading edge (the liquid surface). 

In order to treat the temperature and position dependence in Eqs.. (6) and 

(7), they were evaluated at five different values of axial position, x, up the 

neck of the tank. The x positions were 3.52, 10.6, 17.6, 24.6 and 31.7 in. Each 

axial node in the analysis was assigned one of these x values. The step transi- 

9 

tion from Eq. (6) to Eq. (7) was prescribed by the point at which Gr > 10 . 
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Treatment of conduction and convection by the thermal analyzer was 
straightforward. One aspect of the physical system, however, was not so 
easy to model. This was the heat transfer effect of mass flow of the boiloff 
gas. The measured boiloff rate of the liquid was 0.0508 Ibm/hr. From this 
the heat flux up the vent line due to mass flow can be calculated. 


q(x) = m C p T g (x) 


^0.0508 



Btu \ 
lbm°F' 


T g (x) (°F) 


= 0.0127 T (x) Btu/hr. 

j§ 


( 8 ) 


Neither the 3-D IPM nor any other thermal analyzer available at HREC at the 
time was capable of handling this mass flow option. Therefore, two limiting cases 
were solved to bracket the real problem. At the same time, the Lockheed/ 
Huntsville Computer program was modified to handle Eq. (8). A third case 
was then run on the mofified program to solve the "flowing gas" case. The 
three cases analyzed are: 

1. The venting gas was assumed to be stationary, but its tempera- 
ture constrined only at the liquid surface and the top of the 
next line. 

2. The venting gas was assumed to be stationary, and its temper- 
ature was fixed at the evaporation temperature at all points. 

3. The venting gas was assumed to be flowing out the vent line at 
the rate dictated by boiloff, and its temperature was constrained 
only at the liquid surface and the top of the next line. 

Case 1 

The fallacy in Case 1 is that the heat exchanger effect of the moving 
gas is ignored. The gas nodes reach a steady state temperature in the vent 
line at a value higher than the real case. The heat leakage to the liquid is 
calculated to be higher than the real case. The heat leakage to the liquid is 
itemized in Table 1. The errors are from the uncertainty of HPI thermal 
conductivity. 
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ITEMIZED HEAT LEAKAGE TO LN^ STORED IN MSFC 20-INCH TANK 
WITH 5% GAS FRACTION AND INSULATED WITH 1 INCH OF l/4-MIL 
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The temperature distributions of the gas, the fill line and the vent line 
(tank wall) are shown in Fig. 2. Note that for each case the curve represents 
the fill and the vent line temperatures. Also, for Cases 1 and 2 the gas tem- 
perature is very close to this same temperature. Inspection explains why 
this is the case. The convective heat transfer of the gas is much greater 
than the pure conduction down the gas. Therefore, radial heat transfer is 
much greater than longitudinal with the result that at any longitudinal position 
the temperatures of all components are very nearly the same value. 

Case 2 


The error in Case 2 is obviously that the temperature increase of the 
gas as it ascends the vent annulus is ignored since the gas is held stationary 
and at a constant temperature in the analysis. The heat leakage to the liquid 
is itemized in Table 1. The temperature distributions of the fill line and the 
vent line (tank wall) are shown in Fig. 2. Note that the temperature distri- 
bution is flat in the vent and fill lines in all but the very top of the neck. This 
is due to the constant gas temperature approximation. The real case is ex- 
pected to be somewhere between Case 1 and Case 2. 

Case 3 

This case most closely represents the real situation since the gas 
motion is considered in the heat balance. The resulting heat leak to the 
liquid is itemized in Table 1, and the temperature distribution in the tank 
is shown in Fig. 2. As expected, the distribution for Case 3 fell between the 
two "bracketing" cases. Also, the heat leakage to the liquid for Case 3 fell 
between that of the two others. 

The three calculated values of total heat leakage to the liquid are now 
compared to the test value of 4,38 Btu/hr. 


9 



LMSC/HREC D148760 



(j ) sati^-eaaduiax 


10 


Fig. 2 - Plot of Temperature* vs Axial Distance above the Liquid for the NASA 20-Inch Cryogenic 
Tank Insulated with 1-Inch of Mylar and Foam and 95% Fill with LN_ 
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Case 1: Percent difference = 


7.58 - 4.38 
4.38 


73% 


Case 2: Percent difference = 


4.15 - .4.38 
4.38 


= -5.25% 


5 79 - 4 38 

Case 3: Percent difference = — — -r * — = 32.2% 

4. 3o 


These results indicate that the method used does have a considerable 
effect on the calculated heat rate to the liquid. In this particular comparison, 
although the mass flow case was expected to produce the closest comparison, 
the fixed gas temperature assumption actually produced the closest. Two 
factors should be noted in judging this comparison. First, the accuracy of the 
analysis is based on the accuracy of the assumed thermal conductivity value, 
5x10 ^ Btu/hr-ft-°F. The thermal conductivity of the mylar and foam com- 
posite is dependent on factors such as its preconditioning, its degree of out- 
gassing and its compression. These and other effects can cause an error 
(+20%) in estimating the applied conductivity value of the HPI. The second 
factor is that there is some question concerning the reliability of the test data 
because of changes in the ullage pressure during the test. 

The significant fact is that when a venting tank is analyzed, careful con- 
sideration should be given to the method of modeling the vent gas temperature. 
Different assumptions, although all reasonable, predict substantially different 
heat leaks. 
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2.2 THE NASA/MSFC 105-INCH CRYOGENIC TANK 

The NASA/MSFC 105-in. cryogenic tank was analyzed using the 
Lockheed/Huntsville Computer program. Figure 3 shows a schematic 
of the tank. The liquid level was chosen by setting the gas fraction in the 
top of the tank at approximately 5%. Only the portion of the tank at and above 
the adiabatic separator, including the tank neck, was analyzed. The cryogen 
was LH^ which has a boiling temperature of -423°F. As in the case of the 
20-in. tank, the heat flux to the liquid below the separator was calculated by 
hand. 



where 


k = 5 x 10" 5 Btu/hr-ft-°F 
2 

A s 239.63 ft (calculated from tank dimensions) 
AT = 70 - (-423°F) = 493°F 
AX = 1 in. 


.*. q = 70.8 Btu/hr 

The remaining heat flux to the liquid is calculated in the thermal analysis. 
The analysis was performed for 1 in. of the double aluminized Mylar and red 
polyurethane foam insulation. The thermal conductivity normal to the insula- 
tion was assumed to be 5x10 ^ Btu/hr-ft-°F while that parallel to the insulation 
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was assumed to be 1.55 x10 2 Btu/hr -ft-°F. The analysis includes heat flux 
through the insulation as well as flux down the tank wall, the fill line, the gas 
and the instrumentation wires . Temperature-dependent thermal conductivity 
values for the aluminum tank, the stainless steel and fiberglas fill and vent 
lines, the copper liquid level sensor, the cons tantan lead wires and the hydro- 
gen boiloff gas were used in the analysis. The curve used for the fiberglas is 
shown in Fig. 4. All the heat flow paths in the thermal analysis were the con- 
duction type except the convection paths connecting the gas to the tank wall, the 
fill line and the instrumentation wires . The convection equations used were 
the laminar and turbulent natural convection equations on a vertical wall, Eqs. 
(3) and (4). These were evaluated for -400°F hydrogen gas to obtain: 

Laminar 

h = 0.922(T - 37) 1//4 x _1 / 4 Btu/hr-ft 2 -°F 

Turbulent 

h = 1,17 (T - 37 ) 2 ^ 5 x l f* 


9 

The transition Gr number was fixed at 10 in the analysis. As for the 
20 -in. tank, the two equations were solved at five axial positions which, in 
this case were 4.5, 13.5, 23, 32 and 41 in., respectively. Each axial node in 
the analysis was assigned one of these x values. 


As for the 20-in. tank, three cases were analyzed, two approximations 
and a "flowing gas" case. These three cases were: 


• The venting gas was assumed to be stationary, and its temper- 
ature was constrained only at the liquid surface and at the top 
of the vent line. 

• The venting gas was assumed to be stationary, and its temper- 
ature was fixed at the evaporation temperature at all points. 

• The venting gas was assumed to be flowing out the vent line at 
the rate dictated by boiloff, and its temperature was constrained 
only at the liquid surface and the top of the vent line. 
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Results for all three cases are shown in Figs. 5 and 6 which present 
the temperature distributions for the vent and fill lines, respectively. The 
heat leakage to the liquid is itemized for each case in Table 2. As expected, 
the temperature distribution and heat leakage for Case 3 fell between the two 
"bracketing" cases. 
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Fig. 5 - Plot of Tank Wall and Vent Temperature Distribution for MSFC 105-Inch Cryogenic 
Tank 95% Filled with LH ? and Insulated with 1-Inch Mylar and Foam Insulation 
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Section 3 

CONCLUSIONS AND RECOMMENDATIONS 
3.1 CONCLUSIONS 

Thermal analysis of venting insulated cryogenic tanks requires a proper 
treatment of the convective and mass flow effects of the venting gas. The 
natural convection from the walls of the fill and vent lines, the tank wall and 
the instrumentation wires to the cooler gas results in a heat rate that is 
several times greater than that which would exist with pure conduction between 
the components and the gas. Because this natural convection exists in a 
gravity field, it is essential that it be included in the calculations. 

The heat transport capability of a flowing gas with a temperature gra- 
dient is considerable. This effect must be included in any detailed thermal 
analysis. The actual heat rate associated with mass flow is 

q = m C T. 

If the temperature of the gas flowing into a nodal position is different from 
that flowing out of the nodal position, then there exists a net heat flow due 
to the gas motion coupled with its thermal capacitance. There is a tempera- 
ture gradient in cryogenic vent lines; therefore, there is a net heat flow out 
of the gas node defined by 

% 

Siet ~ ™ ^p ^ T out ^in^ 
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This net heat flow is ignored in a standard thermal analysis using the sta- 
tionary node approximation where the only heat flow in the direction of the 
gas motion is pure conduction in the gas. In the mass flow treatment case 
the vertical heat transport in the gas due to boiloff gas flow out the tank 
vent is approximately one order of magnitude greater than the vertical heat 
flow due to pure conduction alone in the gas. As seen in the comparison 
between cases, the effect of this heat transport on the final temperature dis- 
tribution and the final heat leakage to the liquid is significant. 

It is interesting to note a significant difference between the two tanks. 
For the 20-in. tank (Fig. 2), "the mass flow" (Case 3) temperature distribution 
is nearer the "stationary gas node" (Case 1) temperature distribution, while 
for the 105-in. tank (Figs. 5 and 6), "the mass flow" (Case 3) temperature 
distribution is nearer the "fixed gas temperature" (Case 2) temperature dis- 
tribution. The reason for this is readily apparent. From a gas temperature 
gradient standpoint. Case 1 corresponds to a zero boiloff rate while Case 2 
corresponds to an infinite boiloff rate. Therefore, the greater the boiloff, 
the closer Case 3 will come to Case 2 and the further from Case 1. The 
larger tank, obviously, has the larger boiloff with the result that Case 3 
came closer to Case 2. 


3.2 RECOMMENDATIONS 

The Lockheed/Huntsville Thermal Analyzer Computer program 
is now specifically applicable to cryogenic tank analysis. Laminar and tur- 
bulent, natural and forced convection equations can be used to produce a heat 
transfer coefficient, h, as a function of distance from the leading edge, wall 
temperature, gravitational acceleration and gas type. 

Mass flow equations have been incorporated into the program and are 
properly treating the thermal effect of the moving gas. However, its use 
currently is predicated on a foreknowledge of the mass flow rate. When this 
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parameter is not known from test results, an estimate must made by hand 
calculations. It is recommended that the computer program be further modi- 
fied to incorporate an iterative procedure to arrive at the analytically accu- 
rate value for mass flow rate. Heat rates to the stored liquid would be first 
calculated based on the value of mass flow rate input. A new value could then 
be calculated based on this total heat leakage and the known latent heat of 
vaporization of the cryogen. This iterative procedure will continue until the 
converged value of mass flow rate is reached. 

The Lockheed/Huntsville Thermal Analyzer program has other options 
such as the transient capability. The current program capabilities coupled 
with those recently incorporated or proposed produces a meaningful computer 
tool for application to the storage of cryogenic fluids insulated vented or non- 
vented storage vessels in any g field and under any external thermal condi- 
tions. It is recommended that this computer program be utilized in future 
thermal analysis work dealing with vented cryogenic tanks. 
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